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1. Introduction 

The origin of mass is the most urgent issue of the particle physics today. One of the candidates 
for the theory beyond the Standard Model is the walking technicolor which is the strongly coupled 
gauge theory having a large anomalous dimension y m ~ 1 and approximate scale invariance due to 
the almost non-running (walking) coupling [[T], ^J. The walking behavior is in fact realized in the 
QCD with large number of (massless) flavors Nf which possesses Caswell-Banks-Zaks infrared 
fixed point (IRFP) ^ in the two-loop beta function. The exact IRFP would be washed out by the 
dynamical generation of a quark mass m in the very infrared region \i<m for Nf < NJ, N"' being 
the critical number. However, for Nf very close to NJ, m could be much smaller than the intrinsic 
scale A (3> m), an analogue of Aqcd, beyond which the coupling runs as the asymptotically free 
theory, so that the coupling remains almost walking for the wide infrared region m < |i < A as a 
remnant of the would-be IRFP. The case Nf > NJ is called conformal window, although confor- 
mality is broken in the ultraviolet asymptotically free region beyond A. 

Although the results from the two-loop and ladder approximation of Schwinger-Dyson equa- 
tion analysis [Q] are veiy suggestive, the relevant dynamics is obviously of non-perturbative nature, 
we would need fully non-perturbative studies. Among others the lattice simulations developed in 
the lattice QCD would be the most powerful tool for that purpose. The above two-loop and ladder 
studies suggest that the walking theory if existed would be in between Nf = 8 and Nf = 12. The 
Nf = 8 in particular is interesting from the model-building point of view: The typical technicolor 
model []|] is the so-called one-family model (Farhi-Susskind model) which has a one-family of the 
colored techni-fermions (techni-quarks) and the uncolored one (techni-leptons) corresponding to 
the each family of the SM quarks and leptons. Thus if the Nf = 8 turns out to be a walking theory, 
it would be a great message for the phenomenology to be tested by the on-going LHC. 

Since the pioneering works on the lattice [^J ^] were carried out, a lot of groups have been 
doing lattice studies nowadays. (See Refs. []|] for a review of recent developments.) 

2. Simulation 

2.1 Simulation details 

In our simulation, we use the tree level Symanzik gauge action and the highly improved stag- 
gered quark (HISQ) action without the tadpole improvement and the mass correction in the Naik 
term. It is expected that the flavor symmetry in the staggered fermion and the behavior towards the 
continuum limit are improved by HISQ action. We carried out the simulation by using the standard 
Hybrid Monte-Carlo (HMC) algorithm. We computed the hadron spectrum as the global survey in 
the parameter region and we obtained M n , M p , f % and (yw) as tne basic observable. 

The simulation for Nf = 4 is carried out at j3(= 6/g 2 )= 3.6, 3.7 and 3.8 for various quark 
masses on 12 3 x 16 and 16 3 x 24. We took over 1000 trajectories on the small lattice and about 
600 trajectories on the large lattice in Nf = 4 case. The simulation for Nf = 8 is carried out at 
/$(= 6/g 2 )=3.6, 3.7, 3.8, 3.9 and 4.0 for various quark masses on (12 3 x 32,) 18 3 x 24, 24 3 x 32, 
30 3 x 40 and 36 3 x 48 for various quark masses. We took about 800 trajectories on each size. 
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2.2 Analysis methods 

In this section, we show preliminary results of the chiral perturbation analysis (ChPT) and the 
finite-size hyperscaling analysis (FSHS) in 4- and 8-flavor case. 

If the system is in the chiral symmetry broken phase (#SB), physical quantities, M#, are 
described by the chiral perturbation theory (ChPT); the polynomial behavior. In particular, about 

the pion decay constant f % = F + c\m,f + C2»i^ H . (Here we don't discuss the existence of the 

chiral log.) If F ^ in the above Eq., it is regarded as the ^SB. 

On the other hand, if the system is in the conformal window, Me are described by the finite- 

i 

size hyperscaling relation (FSHS) [|l(]] ; LMh = ^{X) where X = Lm^ y . The 7 in this equation is 
defined as the anomalous mass-dimension. We carry out the hyperscaling analysis with our data of 
M H = {M^,/^,M p } by the following fit function; LM H = c Q + c\X. 
In the following, we analyze Nf = 4 and 8 by these methods. 

3. Spectrum 

3.1 Nf = 4 

In this subsection, we analyze Nf = 4 system by the ChPT and the finite-size hyperscaling 
relation. The result of Nf = 4 is shown in Fig. |], in which the pion mass squared, the decay 
constant and the chiral condensate are plotted on the panel from the left to the right respectively. 
M\ is proportional to rtif. f n and ( \ff\j/) have the non zero value in the chiral limit. Thus, the Nf = 4 
has the property of the ^SB phase and this is regarded as the signal of the chiral broken phase in 
the dynamical case of lattice QCD. Also, if the FSHS test is applied to Nf = 4 which is #SB phase 




Figure 1: In Nf = 4 SU(3) gauge theory on 16 3 x 24 at j3 = 3.7; Left: M\ as functions of rrif, Center: f K 
as functions of m/, Right: (\j/y) as functions of nif. The blue-solid line in the left panel is the linear fit and 
the lines in other panels are the quadratic fit. 

(the ordinary QCD), what happens? The result of this attempt for f n is shown in Fig. ||. From 
these Figs., there is no data alignment in the region < 7 < 2. This is the property of QCD when 
the finite size hyperscaling is applied to. These properties (ChPT and FSHS ) in Nf = 4 may hint 
whether Nf = 8 is ^SB/walking or conformal. 
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Figure 2: Finite size hyperscaling test of f K in Nf = 4 St/ (3) gauge theory; Left: at 7= 0.0, Center: at 
7= 1.0, Right: at 7= 2.0. 



3.2 TV/ = 8 

In this subsection, we analyze Nf = 8 system by the ChPT and the FSHS test. The panels 
in Fig. H are M K , f n and M p at j3 = 3.8 in particular as a function of the quark mass ray, and the 
polynomial fit (quadratic fit) and the power fit are plotted. 

In M % and M p , the plateau appears in mj < 0.06 on small lattice (12 x 32) and in nif < 0.02 
on large lattice (24 3 x 32) at /3 = 3.8. In the corresponding region of the plateau, f n behaves like 
the linear toward the zero. Since these might be the effect of the finite size effect or might be in a 
different vacuum, these data are not included in the following analyses. 

In Fig. |5[ to take the infinite volume limit is difficult. Then we take the data on the largest 
volume at each nif for the fitting. The fit range is 0.0 < rrif < 0.1. The fit result of M 2 . is obtained 
as follows: M 2 = 2.31(2)m / + 12.5(l)m} , (# 2 /dof = 17.9) and M 2 = 5.43(4)m} 197(3) , (^ 2 /dof = 
34.0). For f n , £ 2 (/jr)/dof = 14.7 in the power fit and x 2 {f 7t )/&of = 6.1 in the quadratic fit, then 
f n ->F = 0.0295(3) in the limit m f -> as the quadratic fit result. For M p , % 2 (M p )/&of = 6.5 in 
the power fit and # 2 (M p )/dof = 1.3 in the quadratic fit, then M p = 0.191(8) in the limit nif ->■ 
as the the quadratic fit result. In all cases, since the # 2 /dof in the quadratic fit is better than that in 
the power fit, the chiral limit by the quadratic fit gives the non-zero value of f n and M p . Thus, it 
seems that the Nf = 8 is in the #SB phase. 

Here we discuss the validity of the ChPT fit. We used the expansion parameter defined as 

^= N f(^nr L ) , where F is the pion decay constant in the chiral limit and M K {nif) is the pion 
mass at nif. In our simulation, 3£ ~ 1.2 — 2.5 = 0(1) at the minimum value of M n ~ 0.2 in our 
simulation. Therefore, our result in Nf = 8 is consistent with ChPT. 

Next we consider the chiral condensate by the direct calculation, (*P*P) = Ti[D' h 1 IS q(x,x)], 
and Gell-Mann-Oakes-Renner (GMOR) relation, I = / 2 M 2 /(4m/). In the result of the direct 
calculation, ~ ray /2 and then the chiral limit is very small in the lattice unit. In GMOR 

relation, the solid line is the combination of the quadratic fit results of M 2 and f n . This chiral 
limit is £ ~ 0.0005 in good coincidence with the direct measurement. Thus the chiral limit of the 
chiral condensate is very small in the lattice unit. This means the chiral condensate of the techni- 
quark, (QQ), is tiny value compared with the mass deformation ray . Therefore, in the following, 
we attempt to find the tail (the remnant) of the conformal. 
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Figure 3: ChPT (the quadratic fit) in f n in Nf = 8 SI/ (3) gauge theory; Left: M\ as a function of mj, 
Center: f n , Right: M p . The solid lines are the quadratic fit and the power fit. 



PBP vs. mf GMOR vs. mf , 

Nf=8. P=3.8 Nf=8, solid line <-- fit result of each M n , f n 




Figure 4: Chiral condensate in N f = 8 SI/ (3) gauge theory; Left: (WW) at j3 = 3.8. Right: GMOR relation 
and the line obtained from the quadratic fit results of M\ and f n , at various j3s. 

If the system is in the conformal window, the data are in the good agreement with the FSHS 
having the universal value of y. Because of the tiny value of the chiral condensate, we apply the 
FSHS to Nf = 8 system in order to catch the tail of the conformal if there is. Fig. || is the FSHS 
test of f n for the various y. The data is aligned (collapsing) at around 7=1. 

Therefore, to quantify this alignment, we attempt the linear fit as the leading approximation of 
FSHS. Panels in Fig. [| are the linear fit result of FSHS for M n , f n and M p from the left to the right. 
Although # 2 /dof is not small, the linearity for each observable is reasonably well and it seems that 
the remnant of conformal exists. The result, y{M n ) ^ /(Mp) ^ y{fn) ~ 1-0, indicates the remnant 
of the conformal. This situation is very interesting to construct the walking model. 

4. Discussion and Summary 

We have made simulations of lattice QCD with 4 and 8 flavors by using the HISQ action. 
We obtained the following preliminary result; The Nf = 4 QCD is in good agreement with the 
chiral broken phase and the f % data is not aligned in < y < 2, which is the characteristics of 
QCD applied to the FSHS test. The Nf = 8 is consistent with ChPT (^f -parameter) and is not 
inconsistent with FSHS particularly for the data of the relatively large m/(> F). We extracted the 
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Figure 5: Finite size hyperscaling test of f K in Nf = 8 SU (3) gauge theory; Left: at 7= 0.6, Center: at 
7=1.0, Right: at 7 = 1 .4. 
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Figure 6: Finite size hyperscaling fit by the linear ansatz of f K in Nf = 8 SU (3) gauge theory; Left: M K , 
Center: f n , Right: M p . The open symbol is not included in the fit data. 



7-value from the FSHS test; non-universal 7 and y{fn) ~ 1.0. We show the table for various /Js as 



the very preliminary result (which is revised from that in Ref. [11]). To understand the behavior of 





j3 = 3.6 


j3 = 3.7 


J3 = 3.8 


/3 =3.9 


P =4.0 


ym-M n 


0.64(1) 


0.63(1) 


0.61(1) 


0.56(1) 


0.56(1) 




0.98(2) 


0.99(1) 


0.95(1) 


0.92(1) 


0.91(1) 


7 in M p 


1.02(2) 


0.91(4) 


0.84(3) 


0.79(4) 


0.77(6) 



Table 1: Preliminary. The statistical error only 



Nf = 8, we compare with our simulation of Nf = 12 which is consistent with the conformal [ ]12| ] 



and the Schwinger-Dyson equation analysis on the finite size and mass [13]. According to SD-eq. 
analysis, it is 7~ 1.0 for the near conformal in #SB phase. This might indicate that the system 
with 7 ~ 1.0 is near conformal (walking). Therefore, Nf = 8 would be a good candidate of the 
walking theory. 

We should mention, however, that there are several possible systematic uncertainties not con- 



sidered in this report; As pointed out in Ref. [13], there exists the mass correction in the hyper- 
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scaling relation for the heavy quark region. Then the linear ansatz adopted in Fig. |6] may not be 
sufficient to fit our data. To improve the situation for better understanding, we will accumulate 
more data for various fermion masses and /3s on larger lattices, and carry out detailed analysis 
using those data. 
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